CMB and 21-cm Signals for Dark Matter with a Long-Lived Excited State 



Douglas P. Finkbeiner, 1 Nikhil Padmanabhan, 2 and Neal Weiner 3 

Harvard-Smithsonian Center for Astrophysics, 60 Garden St., Cambridge, MA 02138 
2 Physics Division, Lawrence Berkeley National Laboratory, 1 Cyclotron Rd., Berkeley, CA 94720 
3 Center for Cosmology and Particle Physics, Department of Physics, New York University, New York, NY 10003 

(Dated: August 16, 2008) 

Motivated by the eXciting Dark Matter (XDM) model of Finkbeiner & Weiner, hypothesized 
to explain the 511 keV signal in the center of the Milky Way, we consider the CMB and 21-cm 
signatures of models of dark matter with collisional long-lived excited states. We compute the 
relic excitation fraction from the early universe for a variety of assumptions about the collisional 
de-excitation cross-section and thermal decoupling. Fhe relic excitation fraction can be as high as 
1% for natural regions of parameter space, but could be orders of magnitude smaller. Since the 
lifetime of the excited state is naturally greater than 10 13 s, we discuss the signatures of such relic 
excitation on cosmic microwave background (CMB) and high-z 21-cm observations. Such models 
have potentially richer astrophysical signals than the traditional WIMP annihilations and decays, 
and may have observable consequences for future generations of experiments. 
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I. INTRODUCTION 

To explain the apparent excess of annihilation 
in the Galactic bulge observed by the INTErnational 
Gamma-Ray Astrophysics Laboratory (INTEGRAL) 
2], Finkbeiner & Weiner Q proposed a model of eXcit- 
ing Dark Matter (XDM) in which Weakly Interacting 
Massive Particles (WIMPs) collisionally excite and sub- 
sequently de-excite via e ± emission. This model uses 
the kinetic energy of the WIMP dark matter to create 
pairs, in contrast with light dark matter models in 
which the pairs result from the mass energy of WIMP 
annihilation [e.g.,h|, where the WIMP mass must be less 
than a few MeV 0, @|. Because the XDM WIMP must 
have a weak-scale mass (~ 500 GeV) it retains many of 
the desirable properties of weak-scale WIMPs such as the 
thermal relic freeze-out abundance. 

The lifetime of such an excited state need not be short, 
and indeed, could be of order the age of the universe 
today. This raises the possibility of a long-lived relic 
excited fraction with observable consequences. A simple 
argument shows the large amount of energy potentially 
available from de-excitations - assuming 100% of the DM 
is the XDM WIMP, and the relic excitation fraction is 
Yf , the energy per baryon, p, is 

p = Y f r)M x — = Y f r) m p (1) 

nb Pb 

where r\ is the fraction of the WIMP mass converted to 
kinetic energy by the de-excitation, M x is the WIMP 
mass, and Pdm/ Pb ~ 5. For the fiducial XDM model, we 
take 

r) = (6-2m e )/M x , (2) 

which for mass splitting S w 1.1 — 2 MeV, and M x = 500 
GeV yields r] w 2 x 10~ 7 to 2 x 10~ 6 , or 



This amount of energy, even if inefficiently transferred 
to the gas, could completely ionize the universe many 
times over for Yf — 1/2. For the more realistic case of 
Yf <C 1/2, the consequences depend on when and where 
the energy is deposited, and with what efficiency. 

This paper explores the astrophysical phenomenology 
of XDM WIMP relic excitations. We start by showing 
that for a natural range of cross-sections, the residual 
excited fraction can be high enough (> 10~ 4 ) to have 
measurable consequences. We then explore these conse- 
quences, focusing on the ionization and thermal history 
of the universe, and discuss how observations of the cos- 
mic microwave background and diffuse 21-cm radiation 
might constrain such effects. Our goals here are two- 
fold: to determine whether the specific model of XDM 
proposed to explain the 511 keV excess is constrained by 
other astrophysical probes, and to explore more gener- 
ally the phenomenology of a WIMP with one or more 
excited states. As we shall show, this more general class 
of "XDM" -like models could have a much richer astro- 
physical phenomenology than traditional WIMPs. 

II. KINETIC DECOUPLING AND DECAYS OF 
XDM PARTICLES 

Before addressing the implications of excited states on 
rcionization, we must address two questions within the 
context of the model: how does the kinetic temperature 
of the XDM relate to the photon temperature when de- 
excitation goes out of equilibrium, and what is the life- 
time of the excited state x* ? The former question is 
important for determining the precise value of the relic 
density of while the latter is important for the trans- 
fer of energy from the x* to ionization in the later uni- 
verse. 



p w Y f (l - 10 keV/baryon) 



(3) 
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A. Summary of the XDM model 

The defining feature of the XDM model is that the 
WIMP has an excited state which can be collisionally 
excited, and subsequently decay to e + e~ pairs. The ex- 
cited state could exist due to compositeness of the dark 
matter, or arise from an approximate symmetry of the 
theory. 

For the excited state to be accessible in the Milky Way, 
and relevant for e + e~ production, only a narrow kine- 
matical range must be considered for the mass splitting, 
5. For the decay to the ground state to be energeti- 
cally capable of producing e + e~ pairs, one must have 
S > 1.022 MeV. On the other hand, the kinetic energy 
available for a pair of 500 GeV WIMPs colliding each 
with velocities v ~ 600 km/ s (roughly the escape veloc- 
ity of the galaxy), is 2 MeV, setting an upper bound on 
S. 

To produce a sufficiently high number of positrons to 
explain the INTEGRAL signal, a large cross section is 
required [H, 0] , comparable to the geometric cross section 
set by the characteristic momentum transfer. That is, 
a ~ (M x £) -1 is of the correct size. Such a cross section 
can arise naturally [3J|, but requires the presence of a 
new light scalar </>, with ^ M X S. The x can excite by 
emitting a cj> with amplitude A_ or can scatter elastically 
with amplitude A+. We generally assume A_ ~ A + , but 
this isn't necessary. 

Most of the equilibrium properties relevant to our dis- 
cussion here are ultimately set by the interactions of </>, 
which stays in thermal equilibrium with the standard 
model through its mixing with the Higgs. Thus, the most 
relevant term for the discussion at hand is the </>-Higgs 
coupling 

C D atfrfh. (4) 

When the Higgs acquires a vacuum expectation value 
(vev), this contributes to the mass of the </>. Thus re- 
quiring a tuning better than 1% in parameters yields 
a naturalness upper bound of about a^l0 2 m|/w 2 ~ 
2 x 10 _3 (m 2 i /l GeV 2 ). Assuming a vev, (<f>) ~ m^, 
one finds a mixing angle between the <p and Higgs of 
sin6> w a ra^v/m^ ^ 10 -4 [3!|. Note that the natural 
range of mixing angle is correlated with m^. That is, 
since a< lO 2 ™ 2 ^ 2 , sin (9 < 10 2 m 3 /(100 GeV) 3 . Thus, 
very light (/>'s are naturally more weakly mixed than heav- 
ier (f>'s. 

B. Kinetic Decoupling of XDM 

Although XDM annihilation XX ^ freezes out in 
a fashion similar to usual WIMPs at T ~ M x /20, ki- 
netic decoupling is a somewhat more subtle story. Direct 
elastic scattering xf ~ > xf is both Yukawa and mixing 
suppressed, and is thus inefficient at maintaining kinetic 
equilibrium. 



/ 

/ 

/ 

/ 



s 

\ 

\ 



FIG. 1: Dominant diagram contributing to kinetic equilib- 
rium of X- 
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FIG. 2: Dominant diagram contributing to thermal equilib- 
rium of <f>. 

The dominant process contributing to kinetic equilib- 
rium of x is — > x4> shown in Fig. [1] The scattering 
cross section for this process is 

A 4 

where A is the x ~ X ~ coupling. (We assume A+ ~ A_ 
here for simplicity, although that does not significantly 
change this discussion. Additionally, we do not distin- 
guish between x an d X* &t this temperature T> i.) 
With such a cross section and A ~ 0.1, assuming a ther- 
mal presence of 4>, x w iU remain in kinetic equilibrium 
down to T ~ m^/SO. 

Ultimately, the relevant process for determining the 
decoupling temperature is when <f> decouples from the 
thermal bath. The dominant diagram for this process is 
shown in Fig. [21 Since we are principally interested in the 
lowest possible temperature Td ec , we are interested in the 
situation where the fermion in question is a muon, and 
the cross section (for relativistic \x and non-relativistic <p) 
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is approximately given by 



y* sin 2 9a e 



(6) 



where a cm is the fine structure constant. 

Such a scattering can keep the <f> in equilibrium down to 
below the muon mass for s'm9 = 10" 4 , while for smaller 
mixings the temperature of decoupling is higher (roughly 
1 GeV for sind = 10" 5 , where additional fields, such 
as pions and kaons, are relevant). As such, we limit 
ourselves to the range of Tdec > 100 MeV, although 
one could conceivably stay in equilibrium longer in other 
models. 



Lifetime 



10" 11 sin" 2 9. Thus, for mixing angles sin 8 ~ 10" 4 , 
such particles would decay before nucleosynthesis, and 
in general well before kinetic decoupling of x- The 
most natural region of parameter space (with the lowest 
tuning), 2m^^ra ( j,^ J l GeV, (f> — > /i + /i" and 4> — > tttt 
become available, and the <j) lifetime will range 3 x 
1CT 15 sin" 2 9 < T0/(1 s) < 3 x 10" 18 sin" 2 9. In this range 
of parameters, <f> will have certainly decayed before ki- 
netic decoupling of \- 

Thus, with reasonable values for the mixing parameter, 
we find the scalar is relatively short-lived. However, the 
decay of the excited state x* into the lighter state x will 
naturally occur late in the universe, producing positrons 
and feeding energy into the baryonic fluid. The amount 
of energy will depend directly of the number of relic x* 's 
left over from the Big Bang. 



In addition to the couplings that drive the early ther- 
mal history, the lifetime of both x* and 4> are clearly 
important. The lifetime of x* is crucial, because this 
determines when the energy of the excited state can be 
deposited into the baryonic gas of the early universe. The 
lifetime of <f> is important, as we have ignored its presence 
in the calculations at T « lMeV and we need to see that 
this is justified. 

We begin by considering the lifetime of the x* ■ Both 
excited and unexcited states of the dark matter will come 
into thermal equilibrium in the early universe, with the 
excited state decaying into the lighter state with an ap- 
proximate lifetime 




Consequently, for the parameters under consideration, 
lifetimes in the range of 10 13 s to 10 18 s are quite reason- 
able. 

Although it appears one can make the lifetime much 
longer simply by lowering m^, there is an implicit link 
between and sin 9 because of naturalness. Lower val- 
ues of ffiffi may quickly lead to a highly tuned region of 
parameter space. 

The scalar <f> decays through its mixing with the Higgs, 
and thus has a lifetime 

T0 = sin" 2 0T h (m h = m^). (8) 

That is, the lifetime of 4> is just sin" 2 6 times that of a 
Higgs boson with mass m^. This decay is typically domi- 
nated by a single process. For example, for ^ 1 MeV, 
the 4> decays into 77. A Higgs of this mass has a life- 
time r~3x 10 _4 s Q. Thus, the <j> will have a lifetime 
~ 3 sin" 2 9 x 10" 4 s. Although such extremely light 
4> bosons are potentially interesting, they occupy a very 
tuned region of parameter space. 

For 2m e ^m0^2m At , (f> — > e + e" dominates, and 
the 4> lifetime will range 10~ 9 sin -2 9 < T^/il s) < 3 x 



III. COLLISIONAL XDM FREEZE-OUT 

Excitation and de-excitation of the WIMP proceeds at 
rates Ue and ko, respectively, via reactions of the form 

XX <-> X*X (9) 
XX* ^ X*X* ■ (10) 

For simplicity, we assume that the rates for both channels 
(|9|10p are equal, and neglect double excitations and dou- 
ble de-excitations. We denote the physical (not comov- 
ing) densities of x, x* by n x , n x * and define n = n x +n x * . 
The Boltzmann equation for n x * is then 

— j 1 — \- 3H(t)n x * = —kDn x *n + kEn x n , (11) 

where H(t) is the Hubble constant at time t. Defining 
Y = n x * jn and a dimensionless inverse temperature x = 
S/T, this simplifies to 

dY 

— = -k D n[Y-(l-Y)f(x)e-*] . (12) 

where 

f{x) = ^1 + ^- , (13) 

which is derived for the relation between the excitation 
and de-excitation rates for a suitable approximation of 
the cross section (see Appendix fA"| . 

We assume that the x particles have decoupled from 
the radiation at a much earlier time than is relevant for 
the freeze-out of the excited states, i.e. Td cc 3> 5; the ki- 
netic temperature T and the photon temperature T 7 are 
then related by T 7 w \/T& cc T. The kinetic temperature 
evolves as 

^ = -2H(t)T - ^ [k E n x - k D n x «] . (14) 

The first term describes the adiabatic cooling of the 
WIMPs, while the second is the thermal energy absorbed 
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FIG. 3: The evolution of the excitation fraction with time, for 



Ms 



1 TeV, T dcc = 100 MeV, 5 = 1 MeV and cr„ 



1 (see 

Eq. IA11|) . Solutions are shown for the full coupled Boltz- 
mann equations (solid black). We show the results under 
certain additional approximations as well, in particular with 
temperature coupling turned off (dashed blue), for the approx- 
imation fc_B = kn exp(— 8/T), ignoring the corrections in Ap- 
pendix [X] (red dotted), and for temperature coupling turned 
off and kE = fcc exp(— 8/T) (long-dashed green). 



and injected as the WIMPs excite and de-excite. This im- 
plicitly assumes that the elastic scattering cross-section 
is much greater than the excitation/de-excitation cross- 
sections; this ensures that the kinetic energy gained/lost 
through dc-excitations/excitations is efficiently thermal- 
ized and the x particles maintain a Boltzmann distribu- 
tion. Comparing the RHS of Eq. QT] and Eq. Q31 it is 
possible to simplify Eq. rjj] and obtain 



d(x~ 1 ) + 2H(t) _ _2dY 
dt x 3 dt 



(15) 



Substituting ln(a) = ln(a/ad ec ) for the time variable, 
where adec is the scale factor at kinetic decoupling, yields 
the following coupled Boltzmann equations: 



The evolution of the excited fraction with time (Fig. [3]) 
exhibits the expected features: Y = 1/2 until the tem- 
perature reaches ~ 5, after which time it rapidly falls 
until the Hubble expansion shuts off the de-excitation 
reactions and it asymptotes to its freeze-out value, Yf. 
We also observe that the modifications to the simplest 
formulation - the difference between kE and kr>, and the 
change in the gas temperature due to the changing frac- 
tion of x m the excited state - have comparable effects on 
the freeze-out abundance. Plotting Yf as a function of 
various XDM parameters (Fig. [JJ , we find that a signifi- 
cant residual fraction can survive in some cases, though 
Yf is smaller for parameters favored by the INTEGRAL 
signal [H, 0] , as computed in [|[ . 

While it is convenient to use adec and Xdcc to de- 
fine the initial conditions, there is an alternative which 
yields a useful scaling of Yf with adec and Xdcc- To 
derive this, we start by noting that the combination 
c = («dcc / 'a) a/ 'xdec / x (which compares the de-excitation 
rate to the Hubble expansion) controls the behavior of 
the system. While c> 1, 7 remains at its equilibrium 
value, and the kinetic temperature simply evolves as a -2 . 
This suggests defining the initial condition as the epoch 
when c = Co > 1. This occurs when a 2 = adcc/co, or 
when xo = ^dec^dec/co- If xq <C 1, then Eqs. 1161 and 
[17] will have identical initial starting points if adec^dcc is 
the same, as Cq is just an arbitrary constant. This implies 
that the freeze-out value only depends on the combina- 
tion ckdec^dcc- This will not be true if xq > 1, but in this 
case, the system will remain in equilibrium as the tem- 
perature falls below 5 and the residual fraction will be 
exponentially suppressed to an uninteresting value. The 
combination adec^dec scales with the XDM parameters 
as 



^dcc*^dc 



10 5 5„ 



VlGeV 



1/2 



lOOGeV 



-5/2 



(19) 



Interestingly, it is independent of 8, at least for the as- 
sumption of Tdcc S> <5 made above. This relationship pro- 
vides a better understanding of the scaling of the curves 
in Fig. |4j In spite of the exponential uncertainty in the 
relic excitation, it is natural to have significant and in- 
teresting relic excitation fractions. Thus, we now turn to 
the observable consequences of such a large Yf. 
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where a dcc = k D (a dcc )n(adcc)/H(ad C c)- Using our fidu- 
cial cosmology, we estimate adec 
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IV. OBSERVATIONAL SIGNATURES 

In this section we consider the consequences of WIMP 
de-excitation at early enough times (z > 10) that the col- 
lisional excitation and de-excitation expected for centers 
of halos at late times is unimportant. 

The most obvious observables affected by energy injec- 
tion into the IGM during the "dark ages" (10 < z < 500) 
are the CMB and the high-z 21 cm background. The 
effect of DM annihilations and decays on the CMB has 
previously been discussed [1, [l(J HH ; the case at hand is 
more like DM decay. The energy injection from WIMP 
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FIG. 4: The residual excitation fraction Yj as a function of 
the scattering cross section <r mr , and decoupling tempera- 



ture Td ec , for M x 



500 GeV {solid black) and M x = 1 TeV 
(dashed red). From top to bottom for each mass, the decou- 
pling temperatures are T dcc = 100 MeV, 500 MeV, 1 GeV, 
and 10 GeV. The 511 keV signal favors a mr = 0.1 - 50 for a 
500 GeV WIMP. 



decay, whether from sterile neutrinos [12l ]. superheavy 
dark matter [Tsj], or generica lly (l4j . has been shown to 
affect the ionization history [15], though if the lifetime 
of the decaying particle is longer than the age of the 
universe, perturbations to the CMB power s pec trum are 
small (Til, Ell- Effects on the high-z 21 cm [H [3 and 
structure formation [l9| have also been investigated. As 
we shall see, the XDM WIMP exhibits a much larger 
range of observable signals for natural models, in part 
because a potentially large fraction of the DM can par- 
ticipate, but also because the characteristic energy scales 
are significantly lower and can be absorbed efficiently by 
the gas. 

We begin by estimating the relevant time scales for 
energy deposition, and demonstrate that the e ± pairs 
produced by XDM de-excitations deposit their kinetic 
energy efficiently. We then consider the possible conse- 
quences of this energy deposition and the potential for 
detecting it in current and future experiments. 



A. Energy Deposition Timescales 

The de-excitation of x* deposits energy into the inter- 
galactic medium in the form of non-relativistic electron- 
positron pairs with kinetic energies ~ 100 keV. The dom- 
inant energy loss for such electrons is via collisions, with 



a cross-section for collisional ionization [20j of 
2.23 x 10~ 15 ln(S/13.6) 2 



E 



(20) 



where E is the kinetic energy in eV. The cross-sections 
for excitations and heating are similar. For a 100 keV 
electron, this corresponds to a cross-section of ~ 2 x 
10 -19 cm 2 , or a scattering rate of 



nnPeHV ~ 5 x 10 



-13 



1 



10 



Comparing this to the Hubble time 



H(z) 



10 16 h 



l + z 
10 



-3/2 



(21) 



(22) 



we see that collisional energy deposition is extremely effi- 
cient over the entire redshift range of interest. We assume 
that all the kinetic energy of the electrons is instanta- 
neously partitioned between ionizations, heating and ex- 
citations. 

The above has focused on the deposition of the kinetic 
energy of the e pairs; there is an additional ~ 1 MeV 
available from the rest mass energy of e^. Positrons can 
annihilate to 2 photons at 511 keV, or to 3 photons. The 
37 spectrum of ortho-positronium [2 11 ] is very hard, with 
only 7 x 10 -3 of the power coming out at E 1 < 100 keV 
and 7 x 10~ 6 at £ 7 < 10 keV. At redshift z < 100, 
the universe is nearly transparent to these photons, so 
their energy is effectively lost. At z > 100 the photon 
energy density is high enough that Compton scattering 
happens faster than a Hubble time, so for relevant life- 
times (r x * w 10 13 — 10 14 s) the mass energy of the pair 
must be included [see e.g. [2^, giving rise to a higher 
effective r). 



B. Ionization/Thermal History 

The effects of XDM on the ionization and thermal his- 
tory of the Universe are controlled by two parameters - 
the available energy per baryon e b and the lifetime of the 
excited state r x * . The energy per baryon is determined 
both by the energy splitting S and the residual excitation 
fraction Yt, 



e b = Y f ^(S-2m e c 2 ) 



(23) 



where the n bx are the number densities of baryons and 
XDM particles. To more easily connect to the relic abun- 
dance calculation of the previous section, we fix 5— 2m e c 2 
to 100 keV and n x /n b to 10 -2 below, corresponding to a 
WIMP mass of 500 GeV. The energy per baryon is then 
trivially related to the relic abundance by 



e b « 10 3 Y) eV 



(24) 
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FIG. 5: The ionization fraction x% = n(e~)/n(H) (top), and 
matter temperature (bottom), for various values of the life- 
time, r x * , with Yj held fixed. The baseline scenario, with no 
energy injection from WIMPs, is shown in both panels (thick 
solid line), and T cm 6 is included in the bottom panel (dotted 
line). In all cases we take Yf = 10 _3 e(,. Note that we ignore 
the effects of star formation etc. on the ionization fraction 
and temperature. 

allowing us to express our results in terms of Yf and t x * . 
It is straightforward to relate the results below to cases 
which make different assumptions for the energy splitting 
and number densities. 

We modify the publicly available code RecFast [40j to 
numerically calculate the ionization and thermal histories 
[HH1]. Examples of these for different choices of XDM 
parameters are plotted in Figs. [5] and [6l where we hold 
one parameter fixed while varying the other. 

The effect of varying Yf at constant lifetime is as one 
expects, with an increasing ionization fraction and tem- 
perature as Yf, and therefore e^, increases. The depen- 
dence on the lifetime is more involved. For lifetimes much 
shorter than the age of the Universe at recombination 
(t x » <C 10 13 s), the energy is simply injected into a fully 
ionized medium with no effect. Injecting energy soon 
after recombination r x » ~ 10 s can truncate recombina- 
tion early, resulting in a higher residual ionization level. 
For even longer lifetimes, the injected energy does not 
perturb the baseline recombination, but partially reion- 
izes the Universe. However, unlike the previous case, the 
recombination processes have been shut off by the Hubble 
expansion, and so the ionization fraction monotonically 
increases with time. This behavior qualitatively persists 
for longer lifetimes, but the degree of reionization de- 
creases as t x * becomes a significant fraction of the age of 
the Universe today and there simply has not been enough 
time for the x* to de-excite. 

The gas temperature behaves similarly, except that it 



FIG. 6: Same as Fig. [5] but for various values of Yf with 
t x * = 10 16 s. The fiducial model (Yf = 10" 3 and r x - = 10 16 s) 
is represented by a red solid line in both figures. 

remains thermally locked to the photon temperature un- 
til z ~ 300, much later than recombination. Lifetimes 
shorter than this have little effect on the gas tempera- 
ture. 

The above has concentrated on the homogeneous Uni- 
verse - the formation of collapsed halos could modify this 
in two ways. The first is the x ~ X* collisions could 
de-excite the %*. Such de-excitations only increase the 
kinetic energy of the colliding WIMPs and do not inject 
energy into the IGM. The second effect occurs when virial 
motions inside halos can re-excite the dark matter. This 
only happens for halos with velocities > y/8/M x and 
only becomes significant after the ionization and temper- 
ature have already been considerably modified by star 
formation. 



C. Implications for CMB and 21-cm Observations 

The only two probes of the z > 10 ionization and ther- 
mal history of the Universe are the cosmic microwave 
background (CMB) and the measurements of the 21-cm 
hyperfine splitting in hydrogen. While CMB measure- 
ments are now a mature field, they are less constraining 
because they are only sensitive to the integrated ioniza- 
tion history and do not probe the temperature of the 
IGM. On the other hand, 21-cm measurements probe 
both the temperature and ionization as a function of 
time, but the experimental techniques are less developed 
with the first pathfinder experiments scheduled for the 
near future. We consider both of these in turn below. 
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1. CMB power spectrum 

The polarization anisotropy of the CMB, and its cor- 
relation with the temperature anisotropy, can provide a 
powerful constraint on the ionization history of the uni- 
verse. The CMB polarization is principally induced by 
Thomson scattering [25| . and was first observed on small 
scales by the DASI interferometer [26[ , and on large scales 
by the WMAP [271] . The small-scale polarization is sen- 
sitive to ionization at the epoch of recombination, while 
large-scale measurements probe the epoch of reioniza- 
tion at z « 10. Such data are, in principle, sensitive to 
any perturbation of the ionization history of the universe 
caused by new physics, such as XDM. 

The effect of XDM on the CMB may be conceptu- 
ally separated into two regimes - effects on recombination 
(t x » ~ 10 14 s) and effects on reionization (r x » ~ 10 16 s). 
The dominant effect on recombination may be thought of 
as an increased residual fraction of ionized atoms, which 
broadens the surface of last scattering. The increased 
scattering both washes out the temperature fluctuations, 
and enhances and shifts the polarization power spectra. 
These effects on the CMB were discussed in detail by 
for the case of WIMP annihilation, but the basic physics 
is also relevant here. Fig. plots the temperature and 
polarization power spectra for two examples of XDM pa- 
rameters, with the lifetime chosen to highlight the effects 
on recombination. While the differences in the TT power 
spectrum are degenerate with the slope of the primordial 
power spectrum, these degeneracies are mostly broken 
by the polarization power spectra. Fig. [7] also plots the 
nominal polarization sensitivities of current and future 
CMB measurements. 

The detectability of these changes in the power spec- 
trum will depend on many factors, including degeneracies 
with other cosmological parameters and details of reion- 
ization. Nevertheless, we may estimate constraints on 
the z w 1000 energy injection. The limit on such en- 
ergy injection, marginalizing over the usual cosmological 
parameters, is cdm < 3 x 10 -14 from WMAP, with an 
improved limit of cdm < 10~ 15 eV/s/baryon at z = 1000 
expected from Planck Using Eq. [3] for r\ = 2 x 10~ 7 
we have 

€dm = 10 3 Y f /T x , eV/s/baryon (25) 

yielding constraints of 

Y f < 3 x 10- 4 (t x ./10 13 s) WMAP 
Yf < 10~ 5 (t x »/10 13 s) Planck (26) 

These constraints do not include higher I data from e.g. 
ACBAR [H[ and CBI Using the e a parametrization 
for delayed recombination in (defined by Peebles et al. 
0) Kim & Naselsky find that e a < 0.02 based on WMAP 
and ACBAR [3(J. This constraint also converts to Yf <^ 
2 x 10" 4 (r x ./10 13 s) in agreement with Eq. [13 

For longer lifetimes r x * > 10 15 s, the effect on the 
reionization history (z w 10) could be pronounced. In 




FIG. 7: The effect of XDM on the CMB TT, TE, and EE 
power spectra. The solid (black) line is our fiducial cos- 
mology, with no XDM. The short-dashed (red) line assumes 
XDM with Y f = 1CT 3 , M x = 500GeV, 8 = l.lMeV and 
r x * = 10 14 s, while the long-dashed line assumes Yf = 1CP 2 
with the other parameters the same. For lifetimes much less 
than 10 13 s, the de-excitation of XDM has a negligible effect 
on the CMB, since the Universe is already completely ion- 
ized. For decays much later, the dominant effect is only on the 
largest scales, and therefore hard to disentangle from standard 
reionization. Also plotted are nominal curves for the polariza- 
tion sensitivity in bins of log 10 (^) = 0.05 for the WMAP and 
Planck CMB missions. Uncertainty due to cosmic variance is 
not included in these sensitivity estimates. 

this limit, it is useful to express the information in the 
CMB polarization in the form of a scattering optical 
depth, r, given by 

_ n a T c f 1 da x l (a) 

Ho Ja Ici a 4 y/n m / a s + r> A ' 

where no is the number density of H at z = 0, Xi is 
the ionized fraction of H, a rc f is the scale factor at some 
early reference time, and the Hubble parameter has been 
expressed in terms of Hq, £l m and We are inter- 

ested in the "excess" optical depth due to x* decay, 
At = r(Yf,T x *) — r(Yf — 0), which we compute for 
a range of Yf and r x * (Fig. [SJ. We identify the re- 
gion of parameter space where At is large enough to 
distinguish from the somewhat uncertain standard sce- 
nario (At ~ 0.01) and small enough to be ruled out by 
WMAP (At - 0.1). While the CMB is a good probe of 
relic XDM excitation in some parts of parameter space, 
21-cm experiments have a broader reach. 
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FIG. 8: The excess optical depth for CMB scattering, At 
(see ^IV C The hatched region is accessible to observation, 
with At > 0.1 already ruled out by CMB observations, and 
At < 0.01 difficult to notice in the presence of the reionization 
caused by standard astrophysics. 




FIG. 9: The gas temperature [K] as a function of the freeze- 
out excitation fraction, Yf, and the lifetime, r x *, for two red- 
shifts. 



2. 21 cm observations 

The hyperfine (spin-flip) transition in neutral atomic 
hydrogen provides a source of opacity to the cosmic mi- 
crowave background from the time the gas temperature 
decouples from the CMB (z — 200) until reionization 
(z ~ 10 — 20). In the standard scenario, T$ < T 7 and 
the line appears in absorption. If sufficient energy is in- 




10 100 
1+z 



FIG. 10: 21 cm signals: the fluctuation amplitude for an arbi- 
trary scale, k = 0.04 Mpc -1 (top), and mean (sky-averaged) 
signal (bottom), for various values of the lifetime, r x * , with 
Yf held fixed. The baseline scenario, with no energy injection 
from WIMPs, is shown in both panels (thick solid line), and 
STt = is shown in the bottom panel (dotted line). In all 
cases we take Yf = W~ 3 eb- Line styles and colors are the 
same as Fig. [S]for easy comparison. 

jected by new physics (or any other mechanism) such 
that Ts > Tj, the line will appear in emission. Obser- 
vations of this line have the potential to constrain the 
evolution of the matter power spectrum, and reveal new 
information about the first sources of ionizing radiation 
[TlT [ttI l3l| . Several projects are already underway to 
observe the line at z ~ 7 — 14 (e.g., the Murchison Wide- 
field Array [Ilj], LOFAR Q, and others). The heroic 
observing efforts now underway may be sensitive to the 
expected signal from relic XDM excitation within 5-10 
years. 

Following standard practice, we assume the spin tem- 
perature Ts is well defined in each volume element of 
space, and neglect the subtle variation of Ts with atomic 
velocity [32J needed for high-precision calculations. We 
relate Ts to the ratio of the number densities of ground 
state and excited atoms via 

n go \T S J 

where gi/go = 3 is the degeneracy factor, and T* = 
0.068 K is the temperature corresponding to the energy 
of the transition. 

The mean (sky-averaged) signal 5Tb is given by 

5% = 27(1 - Xi )(l^ (^)mK, (29) 

where Xi is the ionization fraction, standard cosmological 
parameters are assumed, and radial peculiar velocity of 
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FIG. 11: Same as Fig. [10] but for various values of Yf with 
t x * = 10 16 s. The fiducial model (Yf = 10" 3 and r x » = 10 16 s) 
is represented by a red solid line in both figures. 

the gas is neglected [ItJ • Figs. ITUl and [TT1 show 5Tb for 
the same set of parameters as Figs. [5] and [5] 

Computation of the T$ history involves collisional cou- 
pling of Ts to Tk via both H-H and H-e~ collisions 
[]~7L [3^. [3^| . Another important effect is the Wouthuysen- 
Field coupling, in which photons in the Ly a resonance 
region exchange energy with atoms via Doppler shift, and 
also couple to the hyperfme transition via Raman scat- 
tering [HI, [H, H3] • The net result is that the presence of 
Ly a photons more tightly couples Tg to the gas temper- 
ature. We follow standard procedure [UHII to compute 
Ts and the gas temperature, Tk (Fig. [9]). 

Because the galactic synchrotron foreground is orders 
of magnitude brighter than 5Tb, the more relevant signal 
to consider is the signal from fluctuations in the gas den- 
sity. These fluctuations cause the observed signal to vary 
both because the density of the gas (and therefore optical 
depth) varies, and also the ionization state and tempera- 
ture. Following fl7} we combine these effects and present 
the expected amplitude at wavenumber k = 0.04 Mpc -1 
(Figs. [TOl and ITTj) . This scale is arbitrary, but is chosen 
to be roughly the largest scale observable by the current 
generation of arrays. It is clear that at z ~ 10 XDM en- 
ergy injection could provide a substantial enhancement 
over the baseline scenario. However, the first stars or 
other astrophysical sources could also produce a dramatic 
signal, so mapping the history at higher z would be nec- 
essary to unambiguously identify any new physics. 



V. CONCLUSIONS 

Exciting dark matter (XDM) was invented to explain 
the 511 keV signal in the center of the Milky Way. This 
was achieved by converting WIMP kinetic energy into 
excitations, and then into e pairs from subsequent de- 
cays. An unintended feature of XDM is that the lifetime 
of the excited state r* may be long (~ 10 13 — 10 18 s). 
We have computed the expected relic excitation fraction 
Yf left over from the early Universe, and explored a vari- 
ety of observable consequences for various Yf values and 
lifetimes. 

Such features would be generic in a wide class of models 
with excited states, beyond the simple model of [3|. We 
find that for the parameters that explain the 511 keV 
signal (500 GeV mass, cr mr w 0.1 — 50) the expected Yf 
is small. However, for a higher mass particle and weaker 
cross section, very substantial relic densities (up to Yf — 
10~ 2 ) are possible, which may be generic in models with 
multiple excited states. The thermal, ionization and spin 
temperature histories of the universe are sensitive tests of 
the detailed physics of dark matter properties, to which 
collider tests may be insensitive. Upcoming probes of this 
era may show anomalies, giving essential insight into the 
nature of dark matter. 
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APPENDIX A: ESTIMATING (DE-)EXCITATION 
RATES 

The velocity excitation/de-excitation rate coefficients 
(&;£>, &e; scatterings per time per density) are given by 

kn,E(r) = J d 3 Vid 3 V 2 f(vi,r)f(v 2 ,r)(7 D ,E(Vrel)Vrel , 

(Al) 

where /(v, r) is the phase space density of particles with 
velocity v at position r, <TD,E(v re i) is the inelastic scat- 
tering cross-section as a function of the relative velocity 
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v re i = |vi — V2I. Assuming the universe is homogeneous 
and the particles are non-relativistic, the phase-space 
density is given by the Maxwell-Boltzmann distribution, 



/(v,r) = 



m \ 3 / 2 / — ini- 
exp 



2ttT 



'IT 



(A2) 



where m is the mass of the particle, and T is the kinetic 
temperature (in energy units). It is convenient to trans- 
form to center of mass variables, V cm and v re ; which 
decouples the velocity integrals over both particles, 



( m \ 3 / 2 / . 

I / a v re i exp 



AttT) 



AT 



v{Vrel)Vrel 



(A3) 

If, following [3j, we assume the de-excitation cross- 
section, (T mr , is independent of velocity, we obtain for 
the de-excitation rate, 



kn = 4cr„ 



T 



(A4) 



which scales as yT as expected from dimensional argu- 
ments. We approximate the excitation rate (Eq. 2 of Q) 
with 



aVrel = (Tmr V W rei = 4 <5/ TO V^ el > A5 / 171 

v 2 rel < 4(5 /m , (A5) 

where <5 is the energy splitting between \ and \* an d 
we assume the same (velocity-independent) cross-section, 
a mr as for de-excitation. Setting x = v re i \J raj AT yields 

/oo 
dx x 2 e~ x2 ^Jx 2 - 6/T . (A6) 

Changing variables, 

poo 

k E = k D I dy ^e' y ^y~S/T , 



6/T 



(A7) 



which is a known integral (Eq. 3.383 of (38J) giving 

k ° = k °(w) K <(w)^{w) • < A8 » 

where K\ is the modified Bessel function of the second 
kind. Using the fact that zK\[z) ss 1 as z — > 0, we see 
that k E « k D for T > 5. As T < 5, we find 



, /7T<5 /-<5 



(A9) 



This implies that k E decreases slower than the naive 
Boltzmann scaling suggests, although the correction only 
grows as y^S/T. The asymptotic expressions suggest an 
approximation, 



fee « ko\ 1 




exp 



T 



(A10) 



this approximation agrees with Eq. IA8I to within a few 
percent for T > S. 

Finally, it is useful to estimate a numerical value for 
these rate coefficients. As in Q, we assume a mr is deter- 
mined by the momentum transfer, a mr — d mr /8m, where 
Omr is assumed to be independent of S and m. Choosing 
fiducial values, this gives, 



2CT mr GeV" 2 
S 



(lMcY 



with 1 GeV~ 2 = 3.90 x 10^ 28 cm 2 
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